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Abstract 

We analyze four-lepton contact interactions that lead to lepton flavor violating processes, with 
violation of individual family lepton number but total lepton number conserved. In an effective 
Lagrangian framework, the assumption of gauge invariance leads to relations among branching 
ratios and cross sections of lepton flavor violating processes. In this paper, we work out how to 
use these relations systematically. We also study the consequences of loop- induced processes. 

PACS numbers: 
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I. INTRODUCTION 



The Standard Model (SM) of Particle Physics does not allow conversion between lepton 
flavors and thus is not able to accommodate neutrino oscillation phenomena. Therefore, the 
evidence of oscillations (see jl| for a recent review) makes leptons a promising sector where 
to look for clues about new physics. This strongly motivates the quest for other lepton 
flavor violating (LFV) searches; a quest that is already being carried out and that will be 
pushed in the near future by a variety of experiments, including the LHC and the projected 
neutrino factory or the linear collider. 

From a theoretical point of view, LFV originates at high energies and can be conveniently 
described in a model independent way using an effective Lagrangian approach. We will 
construct our effective theory using just SM fields. Under this assumption, the minimal 
extension of the SM that can accommodate neutrino masses consists on adding the following 
dimension five operator to the SM Lagrangian j^: 

C = Csm + ^-(LtHfiLjH) + h.c. (1) 

This operator violates total lepton number, L = ULi = L e + L M + L T , by two units, and 
gives rise to Majorana masses for the neutrinos after the breaking of the electroweak sym- 
metry. This dimension five operator also violates flavor, and induces rates for LFV processes 
suppressed by the scale A. In view of the measured neutrino mass splittings, this scale is 
presumably very high, possibly close to the Grand Unification scale. Therefore, the rates for 
. ^FV processes induced by this operator are probably too small to be observed experimentally 

On the other hand, there could exist violation of just the individual family lepton num- 
bers, ALj ^ 0, while preserving total lepton flavor number, AL = 0. This violation could 
be generated at a lower scale than A in ((TJ. In this paper we will concentrate on contact 
four-lepton interactions, that are unrelated to neutrino masses, and that could generate rates 
for the LFV processes at observable levels. In fact, this turns out to be the case in most 
extensions of the Standard Model, such as supersymmetry. We would like to stress that, 
throughout the paper, by LFV we mean violation of just the family lepton numbers; total 
lepton number is conserved. 

We will assume that the physics underlying the effective theory respects the standard 
SU(3)c x SU{2) L x U(1)y gauge symmetry. (Some of the problems that originate when 
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dealing with non-gauge invariant interactions at low energies have been discussed in Ref . . ) 
This simple and reasonable assumption has a very profound consequence: the same operator 
induces several processes, that are related by gauge invariance. Therefore, one could use the 
power of gauge invariance to translate constraints on well studied processes into constraints 
on more poorly measured processes, or even not measured at all. This rationale has been 
applied in the past to constrain non-standard neutrino interactions 5||, to rule out LFV in- 
fractions as an explanation for the LSND anomaly p or the atmospheric neutrino anomaly 
and to study the prospects to observe LFV in a future muon or electron collider [^|. 

In this paper we will try to be more ambitious, undertaking this analysis for all the 
possible LFV effective interactions, involving only leptons and compatible with the SM 
gauge symmetry. To this end, we will classify all the possible operators that can induce 
LFV processes, and we will list different LFV processes induced by those operators, together 
with their present experimental bounds, if these exist. As explained above, gauge invariance 
relates some of these processes and their corresponding bounds. We will construct tables 
where these relations could be used systematically, in order to translate the bounds on the 
best constrained LFV processes into bounds on the worst constrained ones. 

The paper is organized as follows. In the next section we work out the basis of the LFV 
operators. The tables are introduced in Section III, where of course we explain how to use 
them. In this section we also study the possibility of obtaining bounds from processes where 
our operators enter at one loop. We devote section IV to some comments and to present 
our conclusions. 



II. OPERATOR BASIS 



We are interested in dimension-six operators involving four lepton fields with a 
current x current structure. The restrictions imposed by gauge invariance on these operators 
constructed just with SM fields was studied by Buchmiiller and Wyler in j^. According to 
this reference, four classes of operators can be built involving four lepton fields, namely, 
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A ijM = (La^L^Lk^Li) , (2) 
B m = {Jjfi^a Lj) (L k ^j^a Li) , (3) 

c m = (h^iHi) , (4) 

D m = (l^j){L kllM Li) . (5) 

Here the subindices i,j,k,l refer to the different lepton flavors (e,/i,r). For one of these 
flavors i, the left-handed lepton doublet is denoted by Lj, and the right-handed singlet by 
%. In eq.Q, er = (cr 1 , cr 2 , cr 3 ) are the Pauli matrices. 

There are indeed other structures that can be formed with four lepton fields, like 

(L i j){kL l ). (6) 

As discussed in j^, one can use the Fierz identity 

(1 ± 7 5 U(1 T 1 5 U = ± 7 5 )7 M ]a4(l T 7 5 )7^] cfe , (7) 

so that eq.(jHI) converts trivially to Dkju- 

However, not all the operators appearing in eqs.(J21EI) are independent. In fact, it turns 
out that the set of A and the set of B are linearly dependent. One has to make use of a 
Fierz identity for the Pauli a matrices, 

&ab&cd = ^ad^cb ~ ^ab^cd > (8) 

and a Fierz identity for the Dirac 7 matrices, 

[(1 - 7VU[(1 - 7 5 )7,W = " [(1 " 7 VU(1 - 7 5 )7,U • (9) 
With this, it is easy to prove the relation 

_1 9 \ / A..,, \ 

(10) 

valid for j ^ I. If j = I, then = (no sum over j). 

For our analysis, we have to choose one of the two sets; for convenience we will keep the 
set A and dispose of the set B. Of course, in the context of a specific data set, it may prove 
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convenient to use B instead of A, or even to use all of them keeping in mind the linear 
dependence. 

In each one of the sets, A^m, Cijki, and Dij kh the subindices run over lepton flavor e, /i, r. 
Some operators are lepton family conserving. Even if these are beyond the Standard Model 
too, our aim is to study LFV, so we concentrate in lepton-family number violating operators. 
They may change lepton-flavor numbers by one unit |ALj| = 1 or two units, |ALj| = 2 . 

In order to proceed and select a basis in the operator space, we have to take into account 
that not all possible combinations of flavor indices lead to independent operators. For 
example, trivially A^m = A^uj. The same happens for C, but not for D. In addition, 
complex conjugation does not really lead to a new operator, in the sense that it is obviously 
the real combination g O + g* what appears in the Lagrangian. Finally, in the case 
of the C-type operators, Fierz rearrangements still lead to relations among different flavor 
combinations. 

We have carefully taken into account all these arguments and reached the following conclu- 
sion: in the case of n flavors, there are n 2 (n 2 — 1)/4 LFV A-type operators, n(n 2 — l)(n+2)/8 
LFV C-type operators and n(n — l)(n 2 + n — l)/2 LFV of D-type. This makes a total of 
n(n — l)(7n 2 + 9n — 2)/8 LFV operators. 

It is useful to enumerate the independent operators in the case of two flavors n = 2, that 
for the sake of definition we take as e and /i. We find three A-type operators, two that lead 
to |AL e | = \AL^\ = 1, namely A eeefJi and A emu , and one leading to |AL e | = |AL M | = 2, 
which is A efiefJi . In the C sector, we also find a total of 3 operators that coincide with the 
order we have shown for A. Finally, there are five independent D-operators: D eeefJi , D etMee , 
Dennn, and with |ALj| = 1 and D efJi6tl with |ALj| = 2. Notice our convention that, 

among all possible combinations leading to an equivalent operator, in the first two indices 
we choose the one that has e preceding /i. We will follow the convention of putting the 
first two indices in increasing generation number, namely, e will precede \i and r, and \i will 
precede r. When the two first indices are equal, we choose the two last indices in increasing 
generation number. 

Since there are three flavors in nature, there are in total 18 A-type operators (among 
them 3x2=6 involving actually only two flavors), 15 C-type (with also 3x2=6 involving 
only two flavors), and 33 D-type (with 5x2=10 involving only two flavors). Therefore, the 
basis has a total of 66 LFV operators. 
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The effective Lagrangian describing LFV processes with total lepton number conserved 

is 

£ = ^^JflijkiA-ijki + CijkiCijki + dijkiDij k i + h.c] . (11) 

We have chosen a convenient normalization that will simplify future expressions. The adi- 
mensional coefficients a, c, and d are complex in general. They contain information about 
the high-energy scale A at which the effective Lagrangian arises, 

^f- K F' (12) 

and similarly for c's and aPs. 

III. LFV PROCESSES 

In the previous section we have introduced a minimal basis of operators that define the 
most general LFV effective Lagrangian compatible with the SM gauge symmetry, eqs.(j2HHI). 
For a given process, there are in general several LFV operators in our basis that induce it. 
Conversely, given one operator, it induces in general several LFV processes. This allows 
to relate different LFV reactions. Therefore, the experimental information available on one 
of them could be used to constrain other related reactions. Let us illustrate this with an 
example. 

We consider in detail the exotic /i-decay channel 

fi + -> e + z/ e i/ e , (13) 

that was discussed in in the context of the LSND experiment on u e appearance in a 
beam. Indeed, the LFV decay (fT3*|) could be an alternative explanation of the anomalous 



LSND results 11| without invoking a fourth neutrino. In fact, this is by now past history 
since the present bounds on (|T3|) exclude such an alternative explanation. Independently of 
the LSND experiment, this decay is still of interest to us since it is related by SU(2)i gauge 
invariance to the better measured process \x — > 3e, as discussed by Bergmann and Grossman 
Q], and constitutes a beautiful example of the point that we want to stress in our paper. 

Let us discuss now for this particular example how to derive relations among different 
processes using the effective operators introduced in section II. Clearly, a gauge transforma- 
tion cannot change family lepton number, therefore gauge invariance only relates processes 



with the same number of leptons and anti-leptons of the same generation. For instance, 
the process fi + — ► e + z7 e z/ e , that involves three leptons with \L e \ = 1 and one lepton with 
\L^\ = 1, is related by gauge invariance to fi + — > e + e _ e + that also involves three leptons 
with \L e \ = 1 and one lepton with \L^\ = 1, but not to fi + — > e + z/ e z/ At , since it involves two 
leptons with \L e \ = 1 and two with = 1. This classification is reflected in our tables. 

It can be checked that only the operators A eeefJL and D efJiee contribute to the process 
fi + — > e + z7 e z/ e . Assuming that one of these operators dominates over the other, the branching 
ratio for this process is equal to either |a eee/ J 2 or \d e ^ ee \ 2 , following the normalization chosen 
in eq. (|llj) . (We will come back to the assumption of no fine tuned cancellations later on.) 
Using table I, it is straightforward to find which other processes the operators A eeefl and/or 
D e ^ ee induce . For instance, it can be readily checked that fi + — ► e + e~e + can be induced 
by A eee[li and D e/iee , and also by C eee/1 and D eeefl . The experimental information about these 
processes, when available, is shown in the last column of the tables. In this specific case, 
there is a stringent experimental bound on /i + — > e + e~e + [12] which leads to 



a\ eeil < 0.5 10- 12 



r 2 < 2 irr 12 



die, < 10- 12 • (14) 

For the sake of simplicity in the notation, we will understand here and in what follows 
the modulus squared of a coefficient when just the square is written. Notice that we are 
able to constrain each coefficient separately because of our assumption of barring fine-tuned 
cancellations. As a result, we find that there are other LFV processes severely limited by 
eq. ()14)) . In particular the limit on the branching ratio of the decay yU + — > e + i7 e i/ e coming 
from eq. (jl4|) is O(l0~ 12 ), and is much more restrictive than the direct experimental limit 
obtained in which leads to a 2 eeM < 9 10 4 , d\ < 9 10 4 . All this information can 
be read from table I. Furthermore, from the table we also find that eq. ()14|) leads to strong 
restrictions on e~e + — > e ± n ¥ 1 e~e~ — > e~/i~, etc. We can use our arguments to bound 
processes like v e e — > z/ e /i, which are extremely difficult to observe. 

This procedure can be straightforwardly applied to other processes in the tables. Table I 
contains processes involving the two flavors \x and e, and in the next tables we list processes 
involving r and e (table II), and r and ji (table III). Finally, table IV is devoted to processes 



with the three flavors at work. Each table is subdivided according to the family lepton 
number of the particles participating in the process. 

Let us summarize the use of the tables. First, one has to check the electronic, muonic and 
tauonic lepton number that are involved in the process we want to study and look at the 
corresponding subdivision in the tables. Once we have found the process we are interested 
in, we can read from the table the operators that contribute to that process. Barring 
cancellations, the experimental constraints shown in the table apply to all the operators 
separately, and those constraints on the operators can be translated into constraints on any 
other process in the table. 

Notice in passing that for all the subdivisions, there is at least one process for which 
experimental information is available. However, this is not enough to constrain all the 
relevant operators. Using just the published bounds on different LFV processes, and the 
limits from requiring that the deviation from the SM prediction of the decay rates T(r + — > 
e + u e u T ) and T(r + — > ^u^r) is smaller than the experimental uncertainty, it is possible to 
constrain all the operators, except one operator of the type A, seven of the type C and seven 
of the type D. 

Up to now, we have been considering the effective Lagrangian (jllj) at tree level. We will 
see now that some of the operators mentioned in the last paragraph can be constrained, up 
to an order of magnitude, using loop-induced processes. The idea is that in some of the 
four-lepton vertexes we can take two of the external lines, close them to form a loop, and 
attach an external gauge boson. (Notice that for operators with |ALj| = 2 no loop can 
be closed.) The resulting process is LFV. Take for example /i — > e-y. It can be generated 
through diagrams like the one shown in Fig. ^ 



V e 
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FIG. 1: Diagram for fi — > e'y, (/ = e, fi, r) 

This transition is absent at tree level and thus the loop must be finite (we have checked 
explicitly that the divergence of the diagram cancels). Unfortunately, the precise value of the 
total contribution to the decay cannot be computed, due to the presence of unknown coun- 
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terterms. Nevertheless, one can estimate an order of magnitude in terms of the coefficients 

n 

of the four-lepton interactions [14]. Among all the possible operators in our classification, 
eqs.fllEI), this transition is induced by D e ^ ee , D eee ^ D eflflfl , D eTT ^, and D^ Tre . The 

resulting width is 

r(n -> er/) ~ 10~ 4 a G 2 F m)m\\d^ , (15) 

with d representing any of the couplings of the relevant operators and m t th e mass of the 
fermion running in the loop. The stringent experimental limit on fi — > allows to put 

the following order of magnitude upper bounds. 



d e ee/n d efiee < 10 

d 2 d 2 < 10 
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d^ eT ,d 2 eTTIX < 10- 11 . (16) 

Having fx and e as external lines, we can also attach a Z-boson instead of a photon and 
thus we will have the LFV process Z — > / u ± e =F . The contribution from the loop gives 

T{Z -> /x T e ± ) ~ 10- 6 G 3 F M 7 z \a\ 2 ln 2 ^ (17) 

with a standing for the relevant coupling of any of the three types. The fact that Z is massive 
implies some important differences between the vertexes \xef and iieZ. While \x — > ey is a 
magnetic dipole transition, for the Z the leading contribution to the amplitude is of the type 

prf{F? + FA>f)eZp (18) 

The loop in the effective theory is calculated at q 2 = M| and has the logarithm appearing 
in ([17)1 . Also, while for \x — > ej the external muon and electron must have different quiralities, 
this is no longer the case for (|18p. so there might be now contributions from the three 
types of operators. To be conservative, when finding numerically our bounds, we will set 
ln(A 2 /m 2 ) = ln(k 2 /M 2 z ) equal to 1. 

The list of contributing operators to Z — ► u ± e =F is shown in table W\ The LEP limits 

n 

on this exotic vertex [12J lead to bounds of order 10 on the corresponding modulus square 
coefficients. With the table, we can compare the limits coming from fi — > cf to the cor- 
responding ones from Z — ► fie. When an operator contributes to both (necessarily of the 



D-type) the former are much stronger that the latter (except in the D eree and D eeeT cases). 
However, for the operators contributing to Z —>■ fie but not to, the limits coming from Z 
decays are relevant, although unfortunately not very restrictive. 

The same analysis can be done for \i — r and e — r transitions, mediated by either 7 or 
Z. Again, the relevant information is in table El One can draw conclusions from the table. 
For example one can see that magnetic dipole transitions where a r runs in the loop lead 
to quite stringent bounds. The physical reason stems from the necessary helicity-flip in the 
loop of those transitions. The interested reader can draw other conclusions from the table. 

Let us stress that the limits in table El are just up to an order of magnitude. Despite 
this limitation, we find them very useful, since some operators that could not be constrained 
using tree level processes can be constrained using loop effects. For instance this is the 
case for many operators that involve more than one fermion in the third generation, such 
as A eflTT , C eTTT , D TTefl , etc. At the end of the day, combining results from tree level and 
one loop processes, we find that among the 66 independent operators that contribute to the 
Lagrangian, only very few remain unconstrained: three of type C and two of type D. These 
are C erer , C^ rj[iT , C erflT , D eT£T , and D^ TflT . All of them, except C eTflT , could be constrained 
in the future using e~e~ and jJiTyT linear colliders. 



IV. FINAL COMMENTS AND CONCLUSIONS 



In this section, we would like to comment on two of our assumptions: the absence of 
cancellations and the imposition of gauge invariance. Finally, we present the conclusions. 

When discussing the decay fi + — > e + v e v e we did not allow for fine tuned cancellations. 
Even if unnatural, such cancellations may spoil some conclusions. In the example we pre- 
sented before, the limits in eq. (j!4jl . coming from the \x — > 3e decay, do not hold if there is a 
destructive interference. For instance, one could have two operators, A eeeil and C eeefl , that 
contribute in such a way that 

\o > eeefi ~^Ceee[i\ — 10 > 

while |a eeeA1 | 2 1CT 12 and |c eeeA1 | 2 3> 1CT 12 . This would mean that there would be a 
contribution to the decay ji + — > e + v e v e much larger than our arguments concluded. We find 
this possibility very contrived, although cannot be excluded from first principles. 
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In the future there might be a positive signal in one (or more) of the processes listed 
in the tables. In this case, it will be desirable to push our type of analysis without the 
assumption of fine-tuned cancellations. To do that, we should work with more observables 
than just a branching ratio or a cross-section; we should consider for example polarization 
measurements. Since these involve different combinations of operators, one would be able 
to do a full statistical analysis and bound all the coefficients of the relevant operators. 

As a final comment, we would like to stress that all the relations we have presented in 
the paper are based on the gauge invariance of the operators. Although we find this a very 
reasonable assumption, only the experiments will tell whether the relevant operators respect 
gauge invariance or not. A violation of these relations would be an indication for extra 
effects not considered in the present analysis. For instance, it could happen that for some 
reason dimension-six operators were forbidden, and the lowest dimension operators appear 

n 

at dimension seven or eight |5J. In this case, the relations presented in this paper would not 
hold. 

In conclusion, we have presented a systematic approach to relations between LFV pro- 
cesses related by the SM gauge symmetry. We have restricted ourselves to purely leptonic 
processes. The strategy has been to build the most general effective Lagrangian that pre- 
serves SU{2)l x U(l)y and leads to purely leptonic physics. We have considered operators 
with energy dimension equal to six, which can be rearranged to the shape of current x current. 
We have introduced no new fields neither right-handed neutrinos. The special shape of this 
effective extension has allowed to relate different LFV processes among themselves, and this 
has been applied to the general study of all the 66 operators involving the three flavors. We 
have proposed a systematic framework for studying these relations by means of the use of 
tables where the different LFV processes are listed. We have also discussed the implications 
coming from loop-induced processes. In the tables, one can review the current bounds for 
leptonic LFV searches. 
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TABLE I: LFV processes involving electron and muon flavors. For the coefficients, it should 
be understood modulus squared where just the square is written (for instance, Oeueu stands for 
l a e/ie/J 2 )- Pm-m represents the probability of the transition muonium-antimuonium normalized, 
for the sake of conciseness in the presentation, to Po = 64e 12 [m 3 /Gpm^] 2 = 2.56 1CT 5 with e 
the electric charge. Therefore, Pm-m = Pm-m / Po- Cm the other hand, cr(X + Y — > X' + V) is 
the cross section of the process X + Y — > X' + y' normalized to oq = G 2 F s/3>^, i.e. a(X + Y — > 
X' + Y') = c(X + Y — » X' + Y')/gq. In the tables, near the threshold for heavy lepton production 
in electron- (anti) neutrino scattering, one has to substitute 
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TABLE II: LFV processes involving electron and tau flavors. B{r — > II' I") represents the branching 
ratio of the corresponding rare tau decay, normalized to BR(t~ — > e~P e u T ), i.e. B(t — > ll'l") = 
BR{t — ► III")/ BR(t~ —* e~i? e v T ). For further explanations in the notation, see tabled]. 
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TABLE III: LFV processes involving muon and tau flavors. See explanations in tables HI and ILTI 
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TABLE IV: LFV processes involving electron, muon and tau flavors.. See explanations in tables [I] 
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TABLE V: Operators contributing to LFV processes via one-loop contributions; the resulting 
bounds are nothing but an order of magnitude. 



